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Abstract 

General Relativity (GR), with or without matter fields, admits a natural extension 
to a scale invariant theory that requires a dilaton. Here we show that the recently for- 
mulated massive GR, minimally coupled to matter, possesses a new global symmetry 
related to scaling of the reference coordinates w.r.t. the physical ones. The field en- 
forcing this symmetry, dubbed here quasi-dilaton, coincides with an ordinary dilaton 
if only pure gravity is considered, but differs from it when the matter Lagrangian is 
present. We study: (1) Theoretical consistency of massive GR with the quasi-dilaton; 
(2) Gonsistency with observations for spherically symmetric sources on (nearly) flat 
backgrounds; (3) Cosmological implications of this theory. We find that: (I) The 
theory with the quasi-dilaton is as consistent as massive GR is. (II) The Vainshtein 
mechanism is generically retained, owing to the fact that in the decoupling limit there 
is an enhanced symmetry, which turns the quasi-dilaton into a second galileon, con- 
sistently coupled to a tensor field. (Ill) Unlike in massive GR, there exist flat FRW 
solutions. In particular, we find self-accelerated solutions and discuss their quadratic 
perturbations. These solutions are testable by virtue of the different effective Newton's 
constants that govern the Hubble expansion and structure growth. 



1 Introduction and summary 



An extension of General Relativity (GR) by a mass term, and more general polynomial 
terms, is motivated by the cosmological constant, and dark energy problems. Such an 
extension was thought to be impossible due to the loss of the hamiltonian constraint, 
leading to the existence of a ghost-like degree of freedom, in addition to the conventional 
5 helicity states of a massive graviton [T]; this 6th degree of freedom is referred to as the 
Boulware-Deser (BD) ghost. 

While the BD work established the loss of the hamiltonian constraint for a broad class 
of theories, later it was shown in Ref. [2] that in all massive gravity theories the lapse 
would necessarily enter nonlinearly at the quartic order in fields, suggesting inevitability 
of the loss of the hamiltonian constraint, starting at the quartic order. 

A loophole in these arguments was found in Ref. [3], by pointing out that in a special 
class of theories the hamiltonian constraint may still be present even if the lapse enters 
the hamiltonian nonlinearly. Moreover, Ref. [3] gave an order-by-order Lagrangian free of 
the BD ghost in a particular limit, which prior to that was also thought to be impossible 

Last but not least, in Ref. [Ij the order- by-order Lagrangian of [3] was resummed into a 
diffeomorphism invariant nonlinear theory, which was proposed as a ghost-free candidate 
for massive GR. Moreover, it was shown in [3] that the theory, in the unitary gauge, 
does exhibit the hamiltonian constraint in the quartic order, even though the lapse enters 
nonlinearly beginning from that order. 

The absence of the BD problem in the unitary gauge to all orders was established in 
Refs. O E], providing the proof of ghost-freedom in the full theory. It looks like this 
can also be generalized to a full covariant hamiltonian beyond unitary gauge in which the 
Stiickelberg fields are retained [7], [H]. 

In the Lagrangian formalism on the other hand, the absence of the BD ghost is related 
to a very special structure of the Stiickelberg sector [Hin]. Recently, Mirbabayi [TU] has 
uncovered a number of remarkable features of this sector: (1) He found that there is an 
enhanced symmetry for quadratic fluctuations on an arbitrary background, but in the 
leading order in the strength of the background. (2) He showed the absence of the BD 
ghost in small fluctuations on any background. (3) Last, but not least, Mirbabayi showed 
that the absence of the BD ghost in the decoupling limit is not only necessary, but is also 
sufficient for its absence in the full theory away from the decoupling limit. 

Since its formulation, and starting with Refs. [llj - |19] . massive GR has been used 
to study cosmology, black holes, and other exact or approximate solutions (see Ref. [20] 
for a review of theoretical aspects of massive gravity and Ref. [21] for a review on phe- 
nomenology of general massive gravity theories). One interesting feature of cosmology in 
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massive GR is the absence of spatially flat and closed FRW solutions, while the obtained 
inhomogeneous solutions may still well approximate the observed world [17] . There are 
also self-accelerated solutions for which the metric can be brought to the FRW form at 
an expense of having inhomogeneities in the Stiickelberg fields [T71 [221 [53] . Very interest- 
ingly however, a solution with FRW symmetry does exist for open universe |24t [25] (for 
subtleties on perturbations about these solutions, see [2^ [27] ). 

Massive GR [3] has also been generalized to theories involving more fields. An arguably 
simplest generalization is to introduce a scalar field the value of which sets the graviton 
mass [IT] . One may go further and introduce an additional dynamical tensor field to form 
bi-gravity a la Ref . [28] , but now using the ghost-free theory of Ref . [^ as a basis of the 
construction, as was done in Ref. [29]. Such a bi-gravity, unlike the earlier versions |28] . 
was shown to be free of the BD ghost [29] ; its various classical solutions have been studied 
in |30[ 131] . Furthermore, the tri-gravity generalization was also considered [32]. 

Most interestingly, Hinterbichler and R. Rosen (HR) [33], have recently shown that 
there is a certain overarching order in the tensorial extensions: in D space-time dimensions 
(D > 3) there are at most D gravitons (one of them being massless) that can form 
consistent interacting vertices of the tensor fields. HR proved this by reformulating the 
theory in a vielbein formalism, where they showed that the hamiltonian construction gets 
significantly simplified (see also references on earlier works using the first order formalism 
in Ref. [33]; see Ref. [3l] for further development of the HR construction.). As a bonus, 
by considering compactifications of the Hinterbichler- Rosen module in various dimensions 
one should be able to obtain consistent interacting theories of a finite number of tensors, 
vectors, and scalar^. 

In the present work, we would like to discuss a particular extension of massive GR by 
an additional scalar field. While certain scalar extensions have already been discussed by 
some of us in [T7], here we introduce a special scalar a that gives rise to a certain new 
global symmetry of the Lagrangian. The symmetry transformation involves the scalar 
itself, and the Stiickelberg fields (p"", a = 0,1, 2, 3, that are necessary if one wishes to work 
with a diffeomorphism-invariant action for massive GR. These four fields are scalars w.r.t. 
diffeomorphisms, but do transform under the Poincare group of the internal space of 4'"''s, 
as emphasized by Siegel in [S^. The new global symmetry that we use as a building 
principle for the action involving the scalar a is realized as follows: 

a^a-aMpi, cp" ^ e'^cj)^ , (1) 

^ Owing to the fact that graviton mass scale and the compactification scale need not be related to each 
other, the latter can be taken to be much smaller than 1/m, in which case all Kaluza-Klein modes could 
be decoupled. 
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where a is an arbitrary symmetry transformation parameter. The rest of the fields in the 
Einstein framj§, and the physical coordinates x'^, do not transform. This symmetry fixes 
uniquely, modulo some irrelevant derivative terms, an extension of massive GR by the a 
field; in particular, one consequence of ([T]) is minimal coupling of matter to gravity in 
the Einstein frame (unlike Brans-Dicke theories for instance, which have matter coupled 
minimally to gravity in the Jordan frame). 

To motivate the symmetry ([1]), we recall that massive gravity is built upon a reference 
Minkowski space with the metric 



3. 



where ijab = diag(— 1, 1, 1, 1) is the fiat metric in the internal space of (p°'^s. The fields 
can be regarded as arbitrary non-inertial coordinates of the reference Minkowski spao 
the transformation of (/)"'s in ([1]) is then just a rescaling of the reference space coordinates. 
Thus, o" is a field that enforces invariance when the reference space coordinates get rescaled 
w.r.t. the physical space coordinates, while all the other fields of the theory remain intac10. 

Then, the following natural question arises: what is the relation, if any, of the a field 
to an ordinary dilaton? 

To answer this question, we temporarily transform our action into the Jordan frame 
in Sec. [2j Then, we find that a indeed enters as an ordinary dilaton would in the purely 
gravitational action. Hence, the obtained action is also invariant w.r.t. conventional 
dilatations. In fact, the two global symmetries - dilatations and new global symmetry 
([1]) - are not independent, as far as the pure gravity action is concerned: the latter is a 
linear combination of the former and a global subgroup of diffeomorphisms. However, and 
perhaps not surprisingly, only one of these two global symmetries can be respected once 
the matter field Lagrangian is introduced. We choose this symmetry to be ([1]), using it as 
a guiding principle for constructing our theory with the matter fields. Thus, the obtained 
full theory is not invariant under dilatations, but preserves ([H). 

This is easier to understand by returning back to the Einstein frame: there, the matter 
fields are coupled to the physical metric in a canonical way, with no direct coupling to the 
Stiickelberg fields (j)"". If so, the matter Lagrangian in the Einstein frame cannot be directly 
coupled to the a field in a nonderivative way, since this would violate the symmetry ([T]). 



^We define the Einstein frame in the standard way - the one for which the kinetic term for the graviton 
has the usual Einstein-Hilber form; Jordan frame on the other hand will feature a kinetic mixing between 
the scalar a and the graviton. 

^This does not mean, however, that the theory admits only Minkowski background; nontrivial back- 
ground solutions do certainly exist. 

^Likewise, (^"'s can be regarded as target space coordinates, with the world-volume coordinates being 
a;^, n = 0,1, 2, 3; in this case a enforces that the rescaling of the target space w.r.t. the world-volume be 
a symmetry. 
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This is not the case, however, for an ordinary dilaton, which does couple to the matter 
fields in the Einstein frame without derivatives, thus violating ([T]). Therefore, our a is 
not a dilaton, nevertheless, it is closely related to the latter (is identical to it in pure 
gravity sector), so we will refer to it as "quasi-dilaton" in what follows (likewise, we will 
occassionally refer to ([T|) as "quasi-dilatations" ) . 

Having established these properties, we move on to examine the theoretical consistency, 
and then to deduce the physical consequences of massive gravity with the quasi-dilaton. 

We start out by examining the absence of the BD ghost in the theory. We show that 
the quasi-dilaton turns into the second galileon in the decoupling limit (the first one being 
the helicity-0 mode of the massive graviton [MlEj), and as such it acquires an additional 
enhanced symmetry ~ the field-space galilean invariance: 



where and b are arbitrary constants. As a result, the theory in the decoupling limit is 
free of the BD ghost. We can therefore use Mirbababyi's method to demonstrate that the 
full theory, away from the decoupling limit, is also BD ghost-free. This is done in Sec. 3. 

Consistency of massive gravity with the quasi-dilaton would also be expected if the 
hamiltonian analysis of this theory were to be performed, since the quasi-dilaton does 
not introduce any derivatives in the mass term, and has a canonical kinetic term in the 
Einstein frame. Thus, we expect the constraint uncovered in O E] to be only trivially 
modified in a theory with the quasi-dilaton. 

As we have already mentioned, we introduce couplings to matter in a minimal way so 
that only the tensor field in the Einstein frame couples to matter fields; such a coupling 
does not spoil the consistency of the theory described above. 

As a next step we study phenomenological implications of the model, and in particular, 
whether the theory is capable of recovering the predictions of GR. To this end, we check the 
existence of the Vainshtein mechanism |37^ [38] . Using the bi-galileon Lagrangian obtained 
in the decoupling limit, we show that static sources do exhibit the Vainshtein mechanism 
- the helicity-0 field is suppressed at observable scales. Hence, the model should pass all 
observational tests of GR, at least for (nearly) flat backgrounds. 

Furthermore, we discuss cosmology of the theory. We find that unlike massive GR, 
it admits flat FRW cosmology; in particular, selfaccelerated vacua exist for a broad part 
of the parameter space. We perform a preliminary examination of the fluctuations on 
these de Sitter (dS) solutions and find that, unlike already existing examples in ghost- 
free massive GR, in general both of the scalar degrees of freedom may propagate on such 
backgrounds. 

The paper is organized as follows. In Sec. 2, we introduce the theory of quasi-dilaton 




(3) 
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massive gravity. Section 3 deals with the analysis of the theoretical consistency of the 
model, while the screening of extra scalar forces and cosmology are discussed in Sections 
4 and 5 respectively. 

We conclude the present section by fixing notation and conventions used in the rest of 
the paper. As pointed out already, we use the mostly plus metric convention. The Levi- 
Civita symbol, e^ua/^y is normalized so that £0123 = !> while the symbol with upper indices 
is obtained by using the inverse metric tensor g^'^ to rise the indices, and by multiplying 
the expression by det{g^u), so that £^123 _ _-|^^ Various contractions of rank-2 tensors will 
often be denoted by square brackets as follows: /C^^j = [K], IC^i^lC^^ = [/C^], K.^aK,"^pK,^fj. = 
[/C^], etc. The (ordered) index contractions on the epsilon symbols will be omitted alto- 
gether; for example, e^'^'"^ e'^^pall^u^ai3 = eelin, as well as e^'^'^ Ey^^'Iloip^pcr = £^^£1/1111, 
with an obvious generalization to terms with different number of II's. 



2 Quasi-dilaton massive gravity 

In this section we introduce the ghost-free model of quasi-dilaton massive gravity (QMG). 
The theory is based on massive GR [3l 0] , representing a consistent nonlinear extension of 
the Fierz-Pauli theory [39]. We start out with a brief summary of these theories, followed 
by the construction of QMG, which can be obtained by supplementing the gravitational 
sector of massive GR with the global symmetry ([1]) . 



2.1 Massive GR 

A generic theory of a massive graviton can be written in a manifestly diffeomorphism- 
invariant way via an introduction of four spurious scalar fields, (j)'^{x) [351 HOt I41j . One 
defines a covariant tensor H^^, related to the physical (coupled to matter) metric g^^ as 
follows: 

H^^=g^,^-d^<P'^dA\ab, (4) 

where o, 6 G (0, 1, 2, 3) are internal indices counting the scalars and rjab = diag(— 1, 1, 1, 1) 
is the Minkowski metric in the internal space. In the unitary gauge, the scalars are frozen 
to coincide with the corresponding coordinates (/)"(x) = however it is often more 

helpful not to resort to any particular gauge. 

A covariant action for a consistent massive graviton with mass m - the one that 
propagates (at most) 5 degrees of freedom on a generic background, and in the presence 
of matter, has the following specific form [3]: 
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m 
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+ Smig,ip), (5) 
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where /C is a four-by-four matrix with the elements defined as 



IC^', = 6ii-y^^i^^K¥d^¥v^b, (6) 
and Ui are specific polynomials of the matrix /C 

U2 = 4([/C2] - [ICf) = 2e^a.e'"'-)C^,)C"i, (7a) 
^3 = -[/C]3 + [/C][/C2] - 2[/C3] = e^aje-f'-lC^^lC^^JC^, (7b) 
= + 6[/C2][/C]2 - [/C3][/C] - 3[/C2]2 + 6[/C4] = e^„^pe^^'5^/Ct/C"^/C^,/C'', . (7c) 

Furthermore, Oj's are two arbitrary parameters characterizing a given theory and 5^ is 
the action for matter fields ip, minimally coupled to the metric g^^. The dependence on 
the scalars (p"" comes in the graviton potential involving the /C tensor 0. 

As shown in [3], the defining property of the potential terms ^/—gUi, which is also 
a necessary condition for the absence of the BD ghost in the decoupling limit around 
Minkowski space-time, is that upon substitution 

= 5,. - V,. = 0, nx) = - ird.TT, (8) 

all terms of the type (S^vr)'^, have to collect into a total derivative, thereby rendering the 
potentially dangerous higher-derivative self-interactions of vr nondynamical. This elimi- 
nates the BD ghost from the theory in the DL, leading to a well-posed Cauchy problem 
and a single propagating degree of freedom in vr, describing the helicity-0 graviton. As 
shown by Mirbabayi, the above necessary condition is also a sufficient one for the ab- 
sence of the BD ghost in the full theory [1^. The polynomials Ui{lC) are specifically 
constructed to satisfy this property [3]: upon substitution dS]), the tensor reduces 
to d^diyTT and the antisymmetric structure of the potential makes it manifest that all vr 
self-interactions indeed collect into a total derivative. In four dimensions there are four 
independent such terms. The zeroth order one is just a cosmological constant. The coeffi- 
cient of the quadratic terms normalizes the graviton mass, leaving two free parameters in 
the theory. The Minkowski vacuum corresponds to (j)"" = S^^x^, and the spectrum on flat 
space consists of the five polarizations of a massive graviton, while the potential specifies 
the consistent interactions of those. 

^The square root of a matrix satisfies the property \/3 • \/A — A. In general, there are 16 choices for 
the square root of a 4 x 4 matrix: the correct choice is the one with all four eigenvalues positive, as it is 
implied by the minus sign in the definition of /C^. 
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2.2 Adding the quasi-dilaton 

We would like to promote the purely gravitational sector of the ghost-free massive GR to 
a theory, invariant under the global rescalings of the four scalars cj)"" w.r.t. the physical 
coordinates, (or, to put it in a different way, the dilatations in the internal space). To 
this end, we introduce a canonically normalized field a, and impose the global invariance, 
realized in the Einstein frame as in ([T|). The extended massive gravity action which respects 
this symmetry reads. 



+ j ^'^x^f^Cmig,,,,^}) , 
where we have defined 



(9) 



= 5>i- e^/^^P' ^ g^'^d^rduct^^Vab ■ (10) 

Note that the global symmetry ([T]) constrains the coupling of a to gravity up to deriva- 
tive terms. We will choose this symmetry as a guiding principle for constructing matter 
couplings as well. In particular, we will couple the matter fields to the Einstein-frame 
metric g^^ in the minimal way, without any direct coupling to the Stiickelberg fields (j)"". 
This significantly constrains the interactions of a with matter, allowing only for irrelevant 
derivative interactions, which we will ignore in what follows. 

Note that the global symmetry ([1]) is a linear combination of the global subgroup of 
diffeomorphisms, x^ — )• e~°'x^, and dilatations, which are also a symmetry of the purely 
gravitational sector and are realized in the Einstein frame as 

x^^e'^x^^, a^a-Mpia, ^ ^ ■ (H) 

For constructing the couplings to matter however, we choose to explicitly break the di- 
latation invariance (jlip , retaining ([1]) as the exact global symmetry of the action [§. 

One interpretation of the theory ([9]) is the following. In massive gravity, we introduce 



®In fact, the choice of ([TJ as the guiding principle for couphng the theory to matter can easily be 
motivated from phenomenological considerations. Indeed, we cannot allow any 0(1) coupling of matter 
to cr, if it is to be stabilized/hidden from the experimental tests of gravity. One can easily see this from 
the following reasoning. Massive gravity without the dilaton possesses a built-in property of screening 
extra scalars from observations - the Vainshtein mechanism [37], which originates from the continuity of 
the theory in the m — >■ limit, leading to the agreement of predictions of a massive theory with GR in 
the massless limit. But if a is to couple to matter, the massless limit of ((9} will feature a free dilaton, 
gravitationally coupled to external sources, modifying GR at 0(1). In fact, even without coupling a to 
matter there still is a potential problem of hiding (one combination of) the scalars tt and cr from solar 
system tests. We will however show in Sec. 4 that the Vainshtein mechanism successfully takes care of 
this issue. 
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a fixed reference metric, wliicli is usually chosen to be Minkowski. In QMG, the quasi- 
dilaton appears in the action only through the combination g^" d^(f)'^dp4P [e'^'^^^'^'^'qah) ■, 
so we can think of it as the dynamical conformal mode of a field-dependent reference 
metric. This situation is reminiscent of bigravity theories in which the entire reference 
metric is promoted to a dynamical field and the spectrum involves a massless and a 
massive gravitons. However, this analogy might be somewhat misleading: in bigravity the 
conformal mode would have a wrong sign kinetic term, whereas in QMG the quasi-dilaton 
has a right-sign kinetic term as long as > 0; hence, the spectrum of the theory consists 
of a massive graviton and a massless scalar (see Appendix aJ^. 

One can make a transition to the Jordan frame by performing the following conformal 
transformation 

g^.. = e'"/^'"5M. , (12) 

under which the action becomes 
S.j=j d^x^^sY""'-^ ^1^^'"'^'^^^'^") 



2 



4 



(13) 



Here g^^ denotes the Jordan frame metric with the corresponding notation for the curva- 
ture invariants constructed from it. 



3 Consistency of QMG 

In this section we study the consistency of QMG, which, as we show below, is closely 
related to the consistency of massive GR. The latter theory was constructed to be ghost- 
free in the decoupling limit [31 H], and only later it was shown to be a consistent theory 
at the full non-linear level via Hamiltonian [5], as well as Lagrangian [10] analyses . The 
current understanding is that the decoupling limit predictions in the given class of theories 
are rather powerful and already based on them one can make statements regarding the 
consistency away from the limit jlOj . 

With this in mind, we won't repeat the full canonical treatment of QMG for showing 
the presence of enough number of constraints, required for the absence of ghosts; the 
situation is rather similar to what happens in ghost-free massive GR, since the presence of 

^The term J d^x e^"^'^^^ ^ Aetdi_i4>°-dv<j)a is consistent with all the symmetries and requirements of the 
theory and can be added to the QMG action with an arbitrary coefficient (note that it is already 
present in ([9]) as a part of UiiJC)). However, we will not be exploring this option here (in massive gravity 
without the quasi-dilaton, this term is trivial). 
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the quasi-dilaton field causes minor modifications not affecting the conclusions regarding 
the absence of ghosts in the theory. We will mostly be concerned with the decoupling 
limit analysis, which is best for displaying the physical content. We also comment on the 
robustness of this analysis for making predictions about the unitarity of the full theory 
towards the end of this section. 



3.1 The decoupling limit and Galileons 

As mentioned above, we can gain a better physical intuition about the degrees of freedom 
and nature of their interactions by working in the decoupling limit of QMG. The limit cor- 
responds to zooming to energy scales well above the mass of the graviton, while decoupling 
(at least as much as possible) gravity by sending Mpi — ?• oo. The degrees of freedom we will 
concentrate on are the quasi-dilaton, the helicity-2 polarization of the graviton h^^ and the 
(canonically normalized) helicity-0 graviton vr, defined by (j)^{x) = J^x^ — rj'^^ d^^i: / M-p\inn? . 
We will ignore the helicity-1 polarization A^, since it is guaranteed to carry two degrees 
of freedom (due to an enhanced U (1) symmetry in the decoupling limit) and guaranteed 
not to couple to sources at the linear order. It enters the decoupling limit Lagrangian 
only quadratically, making = a consistent solution to the equations of motion. The 
specific form of the potential y/—gU{lC) guarantees that all dangerous self- interactions of 
the helicity-0 mode, being total derivatives, are rendered non-dynamical and drop out of 
the action, as discussed in the previous section. The decoupling limit is then defined as 
follows: 

Mpi ^ oo, m ^ 0, A3 = (Mpim^) ^''^ = fixed , = fixed . (14) 

In this limit, the Einstein frame Lagrangian ([9]) in terms of the canonically normalized 
fields (with the metric perturbation normalized in the usual way, h^^ — t- /i^^/Mpi), reduces 
to the following expression. 



1 / 3 1\ 1 



+ (7 



1 ^ 1\ 1 _ /I 3 \ 1 04 1 ^^^^ 

-sen + [^-a, - - j ^..nn + (^-a, - -«3 j ^..nnn - y ^.^nnnn 



Mpi ""^^^ 



(15) 
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Here 11 = dndyTi is the rank-2 tensor constructed from the second derivatives of the 
hehcity-0 graviton, T^^ is the matter stress-energy tensor, while 

(-^/i),.. = (nv - d^^d'^K. - d.d'^Kf, + d^d.h - v^,ah + ii^.d'^d^K^ (le) 

is the linearized Einstein tensor on the Minkowski background. 

The decoupling limit theory, despite the presence of higher derivatives in the action, 
is consistent and does not propagate any additional ghost-like degrees of freedom. The 
interactions between the helicity-2 and scalar gravitons come from the potential term in the 
action and, being certain scalar-tensor generalization of Galileons, have been shown to be 
ghost-free in [30] (in the antisymmetric form used above, it is obvious that the presence of 
the Levi-Civita symbol forbids more than two time derivatives on any field in the equations 
of motion). Quite remarkably, the sector of the theory involving the interactions of the 
quasi-dilaton with the scalar graviton is nothing but a bi-Galileon theory, rendering the 
decoupling limit of QMG completely free of ghosts. 

3.2 Comments on the consistency of the full theory 

The unitarity of QMG in the decoupling limit, as the experience with the ghost- free massive 
GR indicates, can be considered as a strong hint of the absence of the BD ghost in the full 
nonlinear theory. The canonical way of obtaining the number of propagating degrees of 
freedom is resorting to the Hamiltonian analysis. Fortunately, the ghost-free massive GR 
has been extensively studied in the Hamiltonian formulation [5] and the presence of just 
enough number of constraints for the propagation of exactly five d.o.f.'s of a massive spin-2 
particle has been shown. It is straightforward to see that the extra quasi-dilaton present 
in the theory can not spoil these constraints, leading to the six propagating degrees of 
freedom in QMG. 

In fact, as shown in jlOj . already the special form of the decoupling limit (jl5p is 
sufficient for proving the absence of the BD ghost in the full theory. At the quadratic 
order, QMG propagates a massive spin-2 state plus a scalar a (see Appendix A) and the 
special Fierz-Pauli structure eliminates the potential BD ghost, which might nevertheless 
show up at the nonlinear order in the pathological interactions of the helicity-0 graviton 
either with itself or the rest of the fields; on generic backgrounds on the other hand, 
if present, the BD ghost will reappear at the quadratic order giving a ghost pole to the 
graviton propagator. The idea is then to study the theory around an arbitrary background 
in a locally inertial frame. Showing the unitarity of the quadratic theory in this setting 
can therefore be considered the proof of the absence of the ghost in the full theory. Let us 
briefly review the argument. The form of the decoupling limit Lagrangian (|15p significantly 
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constrains what the full theory can look like. Based solely on this limit, one can argue 
that the part of the full action involving the "vector" describing the deviation of the 
four scalars from the unitary gauge (f)"" = d^j^x^ — ij'^^A^, has the following schematic form 

£ D h'"' (ai efj_euS + 02 SfiEySS + 03 e^SySSS) 
+ hFF{bi + h2dA + ...) + a{cieeS + C2 eeSS + C3 eeSSS + C4 eeSSSS) ^^^^ 

+ aFF{di + d2dA + ...)- -^FF {ei + €2 dA + . . .) + h^'a"' + ... 

m.n 

Here {ai,bi,Ci,di,ei) are some constants, {m,n) denote arbitrary (semi)positive integers, 
we have defined S^^ = (9(^j4j^), F^^ = and ellipses denote the GR interactions along 

with the (T kinetic term. The key insight of [10], based on the massive GR action without 
the quasi-dilaton, is that on any background characterized by VEVs of the fields appearing 
in ()17p . the specific form of the interactions guarantees that /iqo and /loi perturbations 
represent Lagrange multipliers in the quadratic lagrangian, while one combination of A^i- 
perturbations is non-dynamical, leaving five propagating degrees of freedom. One can see 
from (jl7p that the interactions of a with A^ are just the "scalar versions" of those of 
h^^, characterized by the same antisymmetric structure. Not surprisingly therefore, one 
can straightforwardly see that exactly the same conclusions apply to QMG, leading to (at 
most) six degrees of freedom (the massive spin-2 field and the quasi-dilaton) around an 
arbitrary background - and therefore the absence of the BD ghost in the full theory. 

4 Screening of Extra Forces 

The presence of extra scalar degrees of freedom in QMG is potentially dangerous for phe- 
nomenology, since gravitationally coupled scalars could lead to an unobserved long-range 
fifth force. The same problem is at first sight present in QMG, since o" is a massless field 
which mixes with the helicity-0 graviton and one combination of these two fields will couple 
to matter; however, as we show in this section, this model possesses a built-in mechanism 
for hiding the extra scalars it propagates. The indication of such screening, called the 
Vainshtein mechanism [37] , can already be read off the decoupling- limit lagrangian featur- 
ing the Galileon interactions. Massive GR without the dilaton has been shown to possess 
the Vainshtein screening for a vast part of the parameter space [3l \T7[ [T5] , and so have 
the bi-Galileon theories |42] . Working in the decoupling limit of QMG which represents 
a certain combination of these two theories, we show that the screening mechanism suc- 
cessfully operates in this model as well; around localized sources both the quasi-dilaton 
and the helicity-0 graviton profiles are significantly suppressed within a certain distance 
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the Vainshtein radiuqfl, and the gravitational potential is completely determined by 



the helicity-2 contribution, recovering GR in this region to a very high precisioro. We 
will concentrate on the decoupling limit analysis, since it completely captures all essential 
aspects of the Vainshtein mechanism. This limit is a valid description of physics at dis- 
tance scales A^^ ^ r ^ m^^, where the lower limit marks the regime in which quantum 
corrections become important in the effective theory, while the upper one is due to the 
definition of the DL. Moreover, we will make a particular choice of the parameters, 

as = -4a4, (18) 

for which the interactions between the helicity-2 and helicity-0 gravitons can be com- 
pletely eliminated by a redefinition of the helicity-2 field [3]; as a result, the theory breaks 
into separate helicity-2 and scalar sectors, the former characterizing the GR part of the 
gravitational potential. To this end, we make the following redefinition of the helicity-2 
graviton in the decoupling limit theory (jl5|) . 

V ^ V + ^V + T3^n^.' (19) 

under which the part of the lagrangian involving h^^ reduces to the usual linearized GR 
coupled to the external stress tensor, while the scalar sector lagrangian is given as follows: 

1 / 7_ 7^ \ UJ 

Cs = -^vr een + 2^eenn + -j-^eeUnU - — creeS 

8 \ Ag Ag / 12 

+ -a( een - ^eeUU + 4Seennn - ^eeHnnn ) + -^ttT + ^".^ ttH T^" 
2 V A| A| A§ J Mpi MpiA| 



(20) 

Here we have defined 7± = 1 ± 804; also, S^,^ = d^dyU and T = rf^^T^^ denotes the 
trace of the matter stress-tensor. The above action involves kinetic mixing of the scalars, 
however it will prove to be simpler to study the classical solutions in this non-canonically 
normalized form. 



*Which e.g. for ordinary sources like the Sun is much larger than the galactic scales. 
^We stress again that the helicity-1 mode of the massive graviton is irrelevant in this context since it 
does not couple to the external stress tensor at the linear level. 
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The equations of motion for a and vr respectively are 



- -eeT. + -een - ^eenn + 2^eennn - ^eeimim = , (21) 
6 2 2A3 2A9 > v ; 

4 4A| 2A| 2 A3 

+ 6See5]nn - ^eeSnnn = 3_T . (22) 

Af A| Mpi ^ ^ 

For localized spherically symmetric, static configurations for which the only nonzero com- 
ponent of the stress tensor is Too = M5{r) /r"^ , these equations can be integrated once, 
which, in the spirit of Galileon theories, reduces them to algebraic equations for the first 
radial derivatives of fields 

wA<,-3A^ + 27+A^-4a4A^ = 0, (23) 
^A, + 37_A2 + 72 A3 - 3A, + 47+A,A, - 12q4A2A<, = (^)\ (24) 

where we have used the following notation: 

Ml,m?) ' ""A3r' Air' 

Note that, in this basis, only vr is coupled to sources and the correction to the gravitational 
potential is entirely determined by its profile. One can express in terms of At^ from the 
a - equation, which reduces the system to its final fornj^, 

A,, = - (4a4A^ - 27+A2 + 3A,) , (25) 



407+04 4 48al 5 _ /r,\3 
^^^^-~^^--\-) ■ ^^^^ 

There are a number of observations one can make about the last equation. In the limit 
of large distances, r » r*, the solution is obtained by simply neglecting all nonlinearities, 
Att ~ {r^:/r)'^ and vr has the usual Newtonian profile, modifying the gravitational force at 
0(1) beyond the Vainshtein radius - a manifestation of the famous vDVZ discontinuity 
|43l m]. Moreover, from (|26p one can see that w > 6 is required for the extra scalar force 



^"it is worth to stress once again that we keep here the notation in terms of 7± not to overload the 
expressions; one should however keep in mind that these constants really depend on 0:4. 
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to be attractive in the vDVZ region. Not surprisingly, this coincides with the condition of 
absence of a scalar ghost on Minkowski vacuum, which can be obtained by diagonalizing 
the decoupling limit scalar action (j20p . On the other hand, well within the Vainshtein 
radius r ^ r*, the largest nonlinearity on the l.h.s. (assuming ~ 1) dominates and the 
solution becomes 



while the correction to the gravitational potential is tinj§ 

12/5 



(29) 



In fact, this screening is parametrically larger than in the DGP model [45 1 1461 ITT], leading 
to the impossibility of observing these corrections in the near future. Note that the seeming 
enhancement for a; ^ 1 of the correction to the Newtonian potential in ()27p might be 
misleading; for large values of u the quasi-dilaton becomes decoupled and the predictions 
of the theory should reproduce those of the ghost-free massive GR, for which the choice of 
the parameters given in (jlSp leads to the scalar dynamics governed by quartic Galileon self- 
interactions. This can be directly seen in (j26p . where, for a; — )• oo, the leading nonlinearity 
is precisely that coming from the quartic Galileon, leading to the corresponding screening 
for vr. 

We have therefore shown that within the Vainshtein radius of a point source the grav- 
itational force reduces to an excellent accuracy to the usual Newtonian attraction me- 
diated by the helicity-2 polarizations of the graviton. Vainshtein screening is a robust 
phenomenon and can be expected to work for a generic class of sources. Also, here we 
have considered only a part of the parameter space, but we expect this mechanism to be 
operative for a more general choice of parameters in the theory. We are not concerned 
with a detailed analysis of the solutions at this stage, however what we have shown above 
can already be considered a direct indication of the presence of a built-in mechanism for 
screening extra scalars in QMG. 



^^Although the profile of a is not probed by matter, it is interesting to observe tliat it is also screened 
by tfie Vainsfitein mecfianism, 



^^^^r^/^ (28) 



being significantly suppressed with respect to the Newtonian potential within r 
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5 Cosmology 



Another important property of QMG is the existence of homogeneous and isotropic flat 
FRW solutions. In ghost-free massive GR, the same Hamiltonian constraint that removes 
the dangerous BD ghost from the theory leads to overrestrictive cosmological implications. 
Namely, it constrains the FRW scale factor to be independent of time, leaving Minkowski 
space as the only flat homogeneous and isotropic cosmological solution [T7]. However, in- 
homogeneous/anisotropic solutions can arbitrarily closely approximate the standard FRW 
cosmology for sufficiently small graviton mass, due to the continuity in the m — )• limit 
|17| . The improvement of the situation in QMG is achieved at an expense of the presence 
of the quasi-dilaton in the theory; then, instead of constraining the scale factor to be a 
constant, the extra constraint simply relates its time evolution to that of the quasi-dilaton, 
implying the possibility for nontrivial flat FRW solutions. However, the following natural 
question arises: for cj — )• oo, the quasi-dilaton should decouple from the rest of the fields, 
leaving the dynamics dominated by massive GR; how is then the a; — t- oo limit of the 
homogeneous and isotropic FRW universe of QMG expected to reproduce the inhomoge- 
neous cosmology of the former theory? The obvious clue to this seeming puzzle is that 
there should be a certain kind of obstruction to QMG's flat FRW cosmology for the values 
of Lo, greater than some critial value. We will explicitly see how this works below. 

5.1 Cosmological solutions and self-acceleration 

In the ghost-free massive GR without the quasi-dilaton, as mentioned above, the extra 
constraint that removes the BD ghost also forces the trivial Minkowski space to be the 
only homogeneous and isotropic flat FRW solution. In QMG however, the presence of 
the extra scalar lifts this constraint, giving it a new role of relating the time evolution of 
the scale factor to its own. In order to see how this works, one can concentrate on the 
homogeneous and isotropic field configurations in the Einstein frame theory We start 
from the most general ansatz for flat solutions with this symmetry: 

ds'^ = -N^{t)dt'^ + a'^{t)dx'^ , / = /(t), (f>' = x\ a = a{t), (30) 



15 



and substitute this into 1^ to obtain the minisuperspace actio 



a f da\ uj a? f da\ 
6— ^ +—7^— ( — I 



N \dt J M^^N \dt J 



+ 6m 



+ (l + + 3a4) e^-A^^' (aiV + a^^) - ^ (as + 4a4) e^-/^^^^ (iV + 3a|) 



+ a4e' 



4(t/Mp 



dt 



(31) 



In what follows we assume the matter sector Sm to consist of a perfect fluid with the energy 
density pm and pressure pm- There are four fields in the above Lagrangian, however time 
reparametrization invariance 



git), 



N 



N 



m ' 



guarantees that the a-equation of motion is redundant, being a certain linear combination 
of the rest. Variation w.r.t. / gives the constraint equation 



1 + ^aa + a4 I a^e'^/^^P' 



1 + 2«3 + 3a4 ) a^e^-Z^^Pi 



+ I (as + 4a4) ae^'^^^^P - 046^"/^^^' = k , (32) 



where k is an integration constant. The Friedmann equation can be obtained by varying 
w.r.t. iV, 



g gO-/J\/pi 

(2 + as + a4) - (3 + -as + 3a4) 

4 a 



3 g2(T/Mpi 

(l + -as + 3a4) — 

2 4 



I gSff/Mp 

-(as + 4a4) — 



(33) 



+Pm, 



here and in what follows the dot denotes a derivative with respect to cosmic time, i.e. 
(T = da/{Ndt) and consequently H = d/a = dlna/{Ndt). Notice that eq. (j32p is a fourth 
order algebraic equation for a generic choice of parameters. However a very simple solution 



^^Note that we have retained the lapse N{t) in the action despite the fact that by time reparametrization 
invariance it can be fixed to an arbitrary value as long as f{t) does not equal to one. However, keeping 
it explicitly is quite convenient, since it allows to quickly derive a first-order Friedmann equation for the 
scale factor. 
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exists if /c = 0, in which case the following ansatz can be used 

e-A^pi = ca{t) , ^ = MpiH, (34) 
with c a constant to be determined below. For this ansatz, the Friedmann equation reads 



00 



3 - I ) M^^H^ = 3M|im2 



1 3 

-(as + 4a4)c^ - (1 + ■^'^3 + 3a4)c^ 



+ (3 + -as + 3a4)c - (2 + 03 + 04) 



(35) 



+ Pm : 

while c itself can be obtained from the constraint equation (j32p with k = 
1 + + Q4 j ~ f 1 + ^"3 + 3a4 j c + ^(as + 4a4)c^ - 040^ 



. (36) 



Before analyzing these equations more closely, it is timely to make an observation about 
the Friedmann equation ()35p : its left hand side flips the sign unless oj < 6. This sign 
flip is unacceptable if the theory is to describe real cosmology in the matter dominated 
era in which case the right hand side is strictly positive. Therefore, quite interestingly, 
we are led to conclude that the parameter space leading to self-accelerated vacua capable 
of describing the universe is orthogonal to the parameter space for which the Minkowski 
vacuum is ghost- free. This is also consistent with the expectation that for a; — )■ 00, 
the sigma field decouples and one shouldn't expect to have a homogeneous and isotroic 
solution. 

One obvious solution to ()36p is c = =^ a — t- —00. Substituting this value for the 
quasi-dilaton background into the Friedmann equation, one obtains the conventional FRW 
cosmology with an effective cosmological constant (energy density) 

Ae// = -3M|im2(2 + a3 + a4), (37) 

which is positive for the choice of parameters, such that 

as + Oi < -2. (38) 

A quick inspection of the Einstein frame action ([9]) reveals that at least at the level of the 
background this is a well-defined solution, since (except for the quadratic kinetic term) 
a appears only through e^'^/*'^^' everywhere in the Lagrangian. One should however be 
cautious about the zero expectation value of the latter quantity, since it multiplies the 
entire action involving the perturbations of the auxiliary scalars, pointing towards an 
infinitely strongly coupled vector/scalar graviton sector at the level of perturbations. We 
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will not pursue this solution further in this paper. 
Another solution is given by 

c = 1 , (39) 



which leads to a vanishing cosmological constant in eq. ()35p . One can show from the a- 
equation of motion (see below) that on this solution af asymptotically approaches unity, 
leading to a Lorentz-invariant Minkowski background at late timej^. However, as we have 
noted above, the absence of a scalar ghost requires w > 6 on this background, whereas the 
parameter space leading to a consistent cosmology is exactly orthogonal, w < 6. Therefore, 
the solution for which c = 1 is unacceptable from the consistency point of view and we 
will discard it in what follows. 

The remaining two solutions of the constraint equation are given by the following 
expressions: 

Sees + 8a4 ± J 9a% - 64a4 , 
c.3 = , (40) 



which, when plugged back into the Friedmann equation ()35p . lead to an effective cosmolog- 
ical constant, describing either de Sitter or anti-de Sitter spaces in the absence of matter. 
It is straightforward to analyze these solutions for finding the parameter space with the 
physical de Sitter vacua (i.e. the ones for which both the left and right hand sides of the 
Friedmann equation are positive). The result is 

as/O, 0<a4<^. (41) 

o 

We therefore conclude that for the vast region of the parameter space given by the condi- 
tions (j4ip . there exist de Sitter solutions, capable of describing the late-time acceleration 
of the universe. 

For completeness, we also need to analyze the equation for /. Varying the action w.r.t. 
cr, we obtain 



1 + ^03 + aA (3 + a/) - 2 ( 1 + + 3a4 ) (1 + af)c 



4 J ' V 2 

3 

+ - (ag + 4a4) (1 + 3a/)c2 - 4a4a/c^ = , (42) 



^■^A remark on physical equivalence of different solutions in the theory is in order here. The presence 
of the global symmetry, the quasi-dilatations ([1]), allows us to shufHe any overall normalization between 
the quasi-dilaton a and the auxiliary scalars (^°' . The real invariant, both under the diffeomorphisms and 
the quasi-dilatations, is e^'^^'^'^^g'^" d^4''^di,(j)*'r]ab, leading to physical equivalence of spatially homogeneous 
solutions with the same value for this invarant. It is therefore obvious, that the asymptotics of the solution 
at hand, e'^''*^^' = a, ds^ = —dt^ + a?dx^, (fp — t/a, 0' — is physically equivalent to the usual Minkowski 
vacuum of the theory. 
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which gives / in terms of the scale factor. Using the constraint equation (j36p . this simplifies 
to 

af = l + (3if2 + H) , (43) 

where 

3 \ „ 3 „ \ 3 



K = C 



3(1 + ^as + 04) - 2(1 + ^03 + 3a4 j c + ^ (03 + 4a4) 



(44) 



As already noted above, the equation (|43p shows that in the absence of a cosmological 
constant, the solution approaches the Lorentz-invariant Minkowski background at late 
times, which is unstable for the choice of parameters that can lead to realistic cosmolog- 
ical backgrounds. We do not have such a problem, however, if the late-time solution is 
cosmological constant dominated; this means that the most natural endpoint of the cos- 
mological evolution in the theories at hand is one of the above-described self-accelerated 
de Sitter backgrounds with tiny curvature ~ ?n^. 

5.2 Perturbations 

Here we give some important, preliminary results on perturbations over de Sitter vacua, 
obtained above. To summarize the results on the background evolution, in the (cosmic) 
coordinates used in the previous subsection, the dS solution is (for simplicity, we will set 
Mpi = 1 from now on) 

(j)' = x\ a = \nica), (45) 

a{t) 

where 

C = const = 1 H ;r, (46) 

while c is given in (jlO|) (subject to the conditions ([^T]) for the dS space). 

At the present intermediate stage, it will be more convenient to work in terms of 
conformal time r, transforming to an "almost unitary" gauge in which the background 
metric is g^^ = o?{T)r]^^ and the auxiliary scalars are frozen to their background values 
(f)^ = CT, (j)^ = X*. We define the perturbations of the dynamical fields in this gauge as 
follows, 

9,. = a\%, +h), a = ln(m) + C- (47) 
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The tensor /C, up to quadratic order in perturbations is given by the fohowing expression 



- c (1 + c + -c' + • • ■)^^^ - h^^^xu + hpphp^^xu + . . . : 



(48) 



where all indices are assumed to be raised/lowered with the flat Minkowski metric and 
Ti^y = dP(t)"'dy(t)'^riah = diag{(? , 1, 1, 1). We will also need an expansion to the quadratic 
order of the metric determinant 



1 



1 



1 



-g = l + -h+-h' - -h^'^h^^ + 



(49) 



The Einstein frame action we would like to perturb can be conveniently parametrized by 
separating the pure general relativity sector in the following way 



Se = \ 



+ 



-g [R - 6H^ 



-g 



6H^ 



m 



-g^'^daCfdua 



U20C) + aMi^) + aMiC) 



(50) 



where H denotes the Hubble parameter of a dS solution at hand. The first term describes 
pure GR on dS space (with a CC, consistent with the expansion rate), while the per- 
turbations of the rest of the lagrangian will describe deviation from GR. The quadratic 
perturbations of the second line of (j50p can be written as follows (note that the indices on 
metric perturbations are not raised and prime denotes a derivative w.r.t. conformal time 
t) 



^ J \ (X (X 



2ujH 



hoidiC + (71^00 + 72^m)C 



+ 73^00 + lihoohii + 'j5hoihoi + jehijhij + 77/1: 



2 

a 1 1 



(51) 



With the definitions 



/3o = 4 + 3a3 + 4q4 
/3i = 2 + 3a3 + 604 
/32 = as + 4a4 
/33 = 2 + 03 + q;4 , 
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the coefficients appearing in the two forms of the action are 



71 = 3m^K 

72 = -\m'^c (9cc2/32 - 4(1 + 2c)c/3i + 3(2 + c)/3o) 



4 

; 

4' 



73 = -jH 



2 

m 

74 = 

' 2 

At o 

75 = 

1 + c 

1 if^ 

76 = tt:"^^ (3cc3/32 - 2(1 + 3c)c2/3i + 3(3 + 2c)c/3o - 24/33) - — (^ + 6) 

lb '4 

'^7 = 4"^' (2cc'/3i - 3(1 + c)c/?o + 12/33) + ^(o; + 6) . 
16 8 

The new decoupling limit 

For analyzing the constraints on the parameter space, coming from the requirement of 
the absence of ghosts on the dS vacua of QMG, we resort to the Stiickelberg treatment of 
the perturbations. To this end, it is useful to return to the cosmic coordinates and define 
the metric perturbation in the standard way, (7^,^ = g^J^^ + h^^ , with g^^^ denoting the 
standard background FRW metric, given in (j45p . The (canonically normalized) metric 
perturbation is decomposed as follows, 

K. = K + (52) 

where V is a covariant derivative w.r.t. the background metric and denotes a canon- 
ically normalized vector field, encoding information about the helicity-1 and helicity-0 
modes of the graviton. Furthermore, the spatial part of the vector is decomposed into 
the irreducible representations of the spatial rotation group as follows, 

Ai = Si + ^^b , diSi = 0. (53) 

We will start by examining the perturbation spectrum in the high frequency, decoupling 
limit, defined aa^^l 

m— ^0, H ^ 0, — = u = fixed . (54) 

m 



The GR part of the action is invariant under the substitution (j52p . while the scalar con- 
tribution to the perturbation lagrangian ()5ip . after some rearrangements and using the 



Note that this decoupling hmit difTers from the one considered on Minkowski space above. 
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explicit form for some of the 7-coefficients, is given in the limit (|54p as followa I 

+ 4(76 + 77)(Afe)4 , (55) 



where /O = C + and the following notation for different coefficients has been used 

p = uju'^ + -^, g = K(l--^), r=-^. (56) 
l + c V 1 + cy l + c 

There are a few important observations one can make about the lagrangian for scalar 
perturbations (j55]) . First, we see that w > is required for avoiding a negative kinetic 
term for one of the scalars p. Moreover, another scalar Aq is non-dynamical on the dS 
background! This is a clear manifestation (and a nice consistency check of the calculation) 
of the BD ghost-freedom of the theory - there can only be a single helicity-0 graviton and 
the fluctuations of the quasi-dilaton in the scalar sector propagating on any background. 
Note that none of the scalars are coupled to the external stress tensor in the decoupling 
limit, meaning that in the full theory their coupling to external sources is strongly sup- 
pressed. This is a manifestation of the absence of vDVZ discountinuity. Last but not 
least, integrating out the auxiliary field Aq, one obtains for the kinetic term of the third 
scalar 

5f D (^r-^ydod,bf + ... (57) 



Using the explicit form of the constants (]56p as well as the relation (I46p . one can check 
that, the coefficient in front of the b - kinetic term vanishes! If the coefficient (75 + 77) 
of the last term in (j55p were nonzero, one would conclude that the single scalar degree 
of freedom p remains propagating for w 7^ on any self-accelerating background (at 
least at the quadratic level). Interestingly enough however, one can check by using the 
background equations of motion that this coefficient vanisheJ^. so that this particular 
decoupling limit has nothing to say about the dynamics of b; then this leaves an option 
of two scalar modes potentially propagating on the self-accelerated background. A more 
refined analysis, needed for uncovering this dynamics will be presented elsewherj^. 

As a last step, we remark on one more possible constraint on the parameter space of 



^^In the decoupling limit H54|) . the scale factor can be replaced by unity up to the terms, suppressed in 
this limit. 

^^We thank Claudia de Rham for pointing out that this has to be the case and Andrew ToUey for first 
showing that this indeed is the case. 

^'^The loss of dynamics for various degrees of freedom on non-trivial backgrounds has been already 



encountered a number of times in the context of ghost-free massive gravity models 1181 1251 127] , as well 
as ghosty extensions of the Fierz-Pauli model [48]. But as pointed out above, the present model, at least 
for a certain parameter space, may avoid this problem. 
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a consistent theory. It comes from demanding the absence of ghosts in the vector sector. 
It is straightforward to see (see also j49j . for example) that the relevant condition is 

75 > 0. (58) 

One can check that it does not lead to any further constraints on top of (j4ip . Moreover, 
75 = translates into the following condition on the parameters (consistent with ()4ip ) 

3 

as < —4, 04 = — 1 — -03. (59) 

In this part of the parameter space, the kinetic term for the vector perturbations on the 
dS solution vanishes, leading to no propagating vector modes on the dS vacua of such 
theories. 



6 Conclusions and Future Directions 

We have presented an extension of ghost-free massive GR by a scalar field, dubbed here the 
quasi-dilaton, based on a new global symmetry ([T]) that appears in these theories. The 
symmetry largely determines the form of the action, including the couplings of matter 
to the rest of the fields. The theory enjoys all important properties of massive GR; it 
possesses a Hamiltonian constraint, responsible for the absence of the Boulware-Deser 
ghost; moreover, nonlinear dynamics of the helicity-0 graviton and the quasi-dilaton lead 
to the presence of the Vainshtein mechanism, screening extra forces associated with these 
fields at solar system distances. 

Most interestingly, and unlike massive GR, the theory does admit flat FRW cosmol- 
ogy with potentially interesting phenomenological implications. For a broad class of the 
parameter space, 

0<a4<— , 0<w<6, (60) 
8 

there exist selfaccelerated dS solutions with stable perturbations. Moreover, one combi- 
nation of scalars becomes non-dynamical on these backgrounds (for a special choice of 
parameters, vector modes can also be made non-dynamical). 

Similar to massive GR, the extra scalars do not couple to matter in the decoupling 
limit on self-accelerated backgrounds, implying a strong suppression for these couplings 
in the full theory (and the absence of the vDVZ discontinuity on such vacua). This leads 
to no easily observable fifth forces in the theory at hand. 

Apart from the theoretical constraints, we anticipate a number of phenomenological 
ones as well, bounding different parameters of the theory. One interesting phenomenolog- 
ical implication of QMG for instance is a modified effective Newton's constant, governing 
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the expansion of the universe (different from the effective Newton's constant for the solar 
system processes, or growth of fluctuations); this can be seen e.g. from the Friedmann 
equation (j35p . This would place mild constraints on the parameter oj, coming from the Big 
Bang nucleosynthesis and CMB (essentially from the expansion rate at the corresponding 
times). 

We defer these and other phenomenological aspects of the model to a separate study. 
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Appendix A 

We show in this appendix that at the linear level, QMG propagates the five degrees of 
freedom of a massive graviton plus a massless quasi-dilaton on Minkowski space. Instead 
of adopting the canonical Hamiltonian approach, we will choose a quicker way, showing the 
presence of just enough number of constraints and the corresponding dispersion relations 
describing the desired degrees of freedom. 

At the quadratic order, the Einstein frame Lagrangian Q reduces in the unitary gauge 
to the following expression, 

2 X £2 = -\h--^{£h)^^ - ^ ((V - 2cTv)2 - {h^, - 2cT<)2) - cod'^ad.a. (61) 

It is convenient for the present purposes to work with a (Jordan frame) redefined spin-2 
field, h^^ = h^^ — 2ari^^, in terms of which the Lagrangian rewrites as follows, 

2xC2 = -\h-{8h)^^^ + 2h-'^{d^d,a - vDa) - ^(/ij^ - J?) + (6 - u)d^ad^a. (62) 

The equations of motion for /i^^,^ and a obtained by varying the latter Lagrangian are 
given respectively as follows (we do not distinguish between the upper and lower indices 
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for notational simplicity), 



+A{d^d,a - rj^^Da) = 0, (63) 

{12-2uj)aa = 2{df,d^h^^-ah). (64) 
Taking the divergence of the first equation yields the constraint, 

= 9uh, (65) 

while taking its trace it implies, 

h = 4-5-C7. (66) 

These constraints remove five degrees of freedom from the system, leaving 6 propagating 
modes. Using (|65p and (j66p in the equations of motion, they can be reduced to the 
Klein-Gordon form, 

□o- = 0, (67) 

(°--')(V-4^-)=0> (68) 
indicating that the spectrum consists of five massive and one massless degrees of freedom. 
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